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ABSTRACT 

We present ammonia maps of portions of the W3 and Perseus molecular clouds in order to 
compare gas emission with submillimetre continuum thermal emission which are commonly 
fH \ expected to trace the same mass component in star-forming regions, often under the assump- 



tion of LTE. 

The Perseus and W3 star-forming regions are found to have significantly different physi- 
cal characteristics consistent with their identification as low-mass and high-mass respectively. 
Accounting for the distance of the W3 region does not fully reconcile these differences, sug- 
gesting that there is an underlying difference in the structure of the two regions. 

Peak positions of submillimetre and ammonia emission do not correlate strongly. Also, 
the extent of diffuse emission is only moderately matched between ammonia and thermal 
emission. Source sizes measured from our observations are consistent between regions, al- 
though there is a noticeable difference between the submillimeter source sizes with sources in 
Perseus being significantly smaller than those in W3. 

Fractional abundances of ammonia are determined for our sources which indicate a dip 
in the measured ammonia abundance at the positions of peak submillimetre column density. 
The presence of infrared sources in a large portion of our sample indicates probable internal 
heating and makes any depletion of ammonia in our sources unlikely. 

Virial ratios are determined which show that sources in Perseus are generally not grav- 
itationally bound and that sources in W3 are, although there is considerable scatter in both 
samples. We find that this that external pressure is necessary for cores at small scales to be 
bound while sources and clusters are gravitationally bound on larger scales. 

Our results indicate that assumptions of local thermal equilibrium and/or the coupling of 
the dust and gas phases in star-forming regions may not be as robust as commonly assumed. 
Alternatively, the assumption that ammonia and thermal emission trace the same mass com- 
ponent in these regions may need to be revisited, along with the degree to which the excitation 
conditions within a star-forming region vary. 

Key words: Methods: data analysis - Stars: formation - Stars: protostars - ISM: clouds - 
Radio Lines: ISM. 
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1 INTRODUCTION 

In the study of star formation, the estimation of masses, tempera- 
tures and the internal gas velocity dispersions of protostellar and 
pre-stellar cores is fundamental to understanding the physical state 
of the core and hence to placing constraints on the processes of 
star formation and early stellar evolution. The submillimetre con- 
tinuum emission from cool dust is regarded as the most reliable 
tracer of mass in cores, since its radiative transfer is simple and its 
abundance uncertainties relatively low. However, deriving the dust 
temperatures required for mass calculations from continuum mea- 
surements is not always straightforward, usually because of a lack 
of far-IR data with sufficient spatial resolution and sensitivity to 
match that available in the (sub)millimetre and mid-IR wavebands. 
In addition to this, acquiring kinematic information to determine, 
e.g. the virial state of a core, requires spectroscopy of molecular 
emission lines that trace the denser core gas. 

It is generally assumed that the ammonia inversion transi- 
tions, with critical densities of a few xlO 4 cnT 3 , can be used to 
trace much the same mass component as the submillimetre con- 
tinuum emission and, since the gas and dust should be in ther- 
mal equilibrium at de nsities only a little above this (10 5 cnT 3 e.g. 
iKeto & CaseTTill2008t> . that both gas temperatures and line widths 
from NH3 can be taken as representative of the dust-traced cores, 
allowing the calculation of core masses and virial ratios. 

Several Green Bank Telescope (GBT) studies have been made 
of the NH3 emission from protostellar cores in a number of Galac- 
tic star-forming regions and have been used to derive the impor- 
tant physical parameters of velocity disp ersion and temp erature in 
the dense gas of the core. One of these ( IAllsopdl2012l ) observed 
beam-spaced, 3x3 position grids centred on ~60 cores detected 
in the W3 star-forming giant molecular cloud in the 850 |J.m con- 
tinuum by Moore et al. (2007). A surprising outcome of this study 
is that the cores appear to be poorly defined in NH 3 on scales of 
one arcmin (0.6 pc at the distance of W3), compared to the submil- 
limetre continuum, in which they are compact with radii < 30" . In 
contrast, no significant differences were found between grid cen- 
tre and edge positions in the kinetic temperature, optic al depth, 
colum n density, or velocity dispersion derived from NH 3 JAllsopd 
12012! ;Moore priv comm.). The simplest explanation for this is that 
the NH3 emission is tracing a significantly larger and less dense 
gas comp onent than th e submillimetre continuum. Another major 
study, by iRosolowskv et al. (2008). observed NH 3 towards a se- 
lection of submillimetre-detected cores (among other positions de- 
scribed in the paper) in the Perseus (mainly low-mass) star-forming 
region. Only a single point per core was observed and so no spatial 
information is available in NH3, but a comparison between the two 
studies shows some significant differences between the NH 3 -traced 
cores in the two star-forming regions. For example, the Perseus 
cores have generally lower gas kinetic temperatures, higher NH3 
optical depths and column densities than in the high-mass star- 
forming regions of W3. Further, if we use NH 3 kinetic tempera- 
tures and velocity dispersions to calculate the ratio of the mass to 
the virial mass (the virial ratio) of cores in the two regions, we find 
that the median val ue in W3 is clo s e to u nity. This result is similar 
to that obtained bv lDunham et alj | |2010|) for the Gem OBI high- 
mass star-forming region, but the median virial mass in the Perseus 
sam ple, when calcul ated in the same way, is a factor of ~10 lower 
(see lKirk et al .120071) . Since we might expect virial ratios to cluster 
arou nd unity on prob ability grounds, the Perseus result is surpris- 
ing. l lKirk et al.ll2007l) find that virial ratios in Perseus are typically 
closer to unity when pressure confinement effects are taken into 



account. However, the disparity in virial ratios calculated using a 
common method may indicate that the NH 3 is not tracing the same 
mass component as the submillimetre continuum data used to ob- 
tain the core masses and radii. Since larger gas volumes should 
generally have larger velocity dispersions, we might be able to ex- 
plain this puzzle, as well as the lack of definition in the W3 cores, 
by a discrepancy in the densities traced by each of NH 3 and the sub- 
millimetre continuum and a different density gradient in the cores 
of the two regions, linked to differences in low- and high-mass star 
formation. 



2 OBSERVATIONS 
2.1 Observational set up 

Fully-sampled maps in the NH 3 (1,1) and (2,2) rotational inver- 
sion transition lines were made using the newly-commissioned K- 
band focal plane array (KFPA) in shared-risk time on the 100-metre 
Green Bank Telescope (GBT), operated by the National Radio As- 
tronomy ObservatorjQ. 

Observations were made over seven sessions from the 16 lh of 
December 20 1 to the 28 th of March 20 1 1 . Sky subtraction was at- 
tained through off-source measurements and flux calibration on the 
scale was achieved through switching observations of a noise 
diode. Comparison with previous observations using the dual-feed 
K-band receiver indicate that absolute flux calibration is accurate 
to <20%. 

Atmospheric opacity values (at zenith) were determined from 
local weather models, collated in the archive maintained by 
R.Maddalena. A spectral bandpass of 50 MHz was used, incorpo- 
rating both the NH 3 (1,1) and (2,2) rotational inversion transitions 
at ~23.69 GHz and ~23.72 GHz respectively. 

Each map was completed using the 'Daisy' pattern, consisting 
of a petal-shaped scan trajectory with continuous sampling. The 
oscillation period of the daisy scan is calculated for each map so 
that full sampling is achieved. The integration time per map pixel 
is radially dependent, resulting in reduced signal-to-noise ratios to- 
wards the map edges. This is compensated for in the reduction pro- 
cess via weighting by effective integration time per pixel. The fidu- 
cial radius of each of our maps is 3.5', with some coverage (incom- 
pletely sampled) beyond this radius due to the size of the receiver 
array. Where necessary, some maps were mosaicked together so 
that complete, fully-sampled coverage of all desired regions could 
be obtained. Sampling rates were set at 1 or 2 seconds per sample, 
dependent upon the source and observing conditions, where the os- 
cillation period of the daisy scan was defined with four samples 
per beam. Typical (median) pointing offsets were ~4" and weather 
conditions were stable during all seven observing sessions. On two 
occasions (the 13 th of lanuary and the 7 th of March 201 1) snow in 
the GBT dish was a concern with regard to observing efficiency. 
Through calibrator observations and comparison with previous ob- 
servations gain values were adjusted where appropriate and we 
have achieved an absolute flux calibration accuracy of better than 
20%. 

Maps were centred on individu al submillimetre continuum 
sources selected from the catalogues of lHatchell et al.l J2005ll2007t) 
for Perseus and lMoore et al. I J2007h for W3. Targets were chosen 

1 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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partly at random but also selected to include both clustered regions, 
so that single maps cover several neighbouring objects, and isolated 
sources. The degree to which the resulting samples are representa- 
tive of the populations of both regions in examined in Section [3~3l 
Thirteen maps were completed in Perseus and 15 in the W3 GMC, 
with some overlaps. The pointing centres, observation dates and 
total integration times of these maps are listed in TableQ] 



2.2 Data Reduction and LTE Analysis 

Individual feeds, i.e. time series data from each of the 7 receivers, 
were reduced by running data through the GBT pipeline which cal- 
ibrates the data and reconstructs spatial image cubes using refer- 
ence scans and performing sky subtraction as appropriate The 
output cubes were produced with an angular pixel size of 6" (in 
comparison to the GBT beam size of 30")- The GBT pipeline uses 
Parseltongu^l scripts to control AIPS tasks which produce images 
of the combined spectral data and removes a second order spectral 
baseline. The mosaicking and weighting of overlapping mapping 
observations was performed as part of the GBT pipeline reductions. 

The determination of optical depth, rotational temperature, ki- 
netic temperature and excitation temperature from ammonia (1,1) 
an d (2,2) spectra , assum ing LTE, has been described in detail 
by iHo & Townesl Jl983 ) and applied to o b servational data by 
numerous authors (e.g. Morgan et alj l20ld : iFriesen et alj 1 20091 ; 
iRosolowskv etalj|2008l). The exact pr ocesses used here are de- 
scribed in detail bv lMorgan et al. I d2012l) . We fit the hyperfine struc- 
ture of the ammonia (1,1) and (2,2) lines using our own developed 
ID10 routines in conju nction with the IDL astronomy library rou- 
tines jLandsmaiJll993l) . The results of the fitting allowed us to pro- 
duce maps of kinetic and excitation temperature, linewidth, opti- 
cal depth and column density. It should be noted here that our de- 
duction of ammonia column density makes use of the assumption 
of local thermodynamic equilibrium (LTE), i.e. that the excitation 
temperature (r cx ) of a source is equivalent to its kinetic tempera- 
ture (7\). The fact that derived values of T CK are not approximately 
equal to values of 7V reflects the fact that T CK is highly dependent 
upon the filling factor of the source in question as it is derived from 
a single transition.. The filling factor of each source has bee n cal- 
culated from the assumption of LTE (see lMorgan et al.ll2010h - 

A portion of the W3 GMC has previously bee n observed and 
mappe d in the (1,1) and (2,2) lines of ammonia bv lTieftrunk et alj 
i 19981) . In addition. IRosolowskv e t al. (2008) performed a single- 
pointing survey of ammonia towards selected positions in the 
Perseus molecular cloud. A comparison of the derived physical pa- 
rameters found by those authors, at positions for which we have 
matching data, shows consistency between the measurements at a 
level reflecting measurement errors of 10-20% (see Sectioni}. 



here. Individual NH3 (1,1) integrated intensity images are presented 
in Figure [3] Noise values vary slightly between the Perseus and 
W3 regions, reflecting the lower integration times for some of our 
Perseus data. Individual r.m.s. values are presented with each map, 
but the Perseus data have r.m.s. off-source integrated intensity val- 
ues in the range 0.36 - 0.86 K km/s with a mean of 0.49 K km/s. 
For W3 the range is 0.18 - 0.47 K km/s with a mean of 0.28 K 
km/s. The radial variation in noise due to the inconsistent sampling 
of the 'daisy' pattern observing mode is obvious in Perseus map 1 1, 
the reader is reminded that only the inner portion of each map, to a 
radius of 3.5', is completely sampled. 

The ammonia maps sample the majority of the Perseus and 
W3 regions associated with submillimetre emission. Although the 
sampling can be seen to be less than complete, the brightest, most 
extended sources in the submillimetre are covered by our obser- 
vations. Our observations of Perseus cover at least some portion 
of each of the well-known regions identified as B5, IC348, Bl, 
NGC13 33 andL1455/L144 8. The High Density Layer (HDL) iden- 
tified bv lLada et all dl978l) on the eastern edge of the W3 GMC is 
covered approximately equally with the more diffuse material in 
the southwest portion of the cloud. The W3 ammonia maps cover 
1178 arcmin 2 split approximately equally (51-49%) between the 
HDL and the more diffuse region. Despite the equivalence in area 
between the two regions, more sources are found in the HDL. Our 
maps cover 1 92 of the submillimetre sources from the catalogue of 
iMoore et all d2007l) . of these, 127 are within the HDL and 65 are to 
the southwest of the W3 region. 

The integrated intensity maps presented in Figure [3] show a 
mixture of condensed cores (e.g. Perseus Map 08) and extended 
structure (e.g. W3 Map 13) with multiple examples of apparent fil- 
amentary morphology (e.g. Perseus Map 06, W3 Map 14). A large 
proportion show multiple cores contained within apparently con- 
nected extended structures (e.g. Perseus Map 12). 



3 RESULTS 

Figure s [T] and [2] show t he sub millimetre maps of lHatchell et alj 
d2005h and lMoore etafl d2007l) covering the Perseus and W3 star- 
forming regions respectively. Overlaid on the images are circles 
depicting the coverage of the GBT-KFPA ammonia maps presented 



2 https://safe.nrao.edu/wiki/pub/Kbandfpa/KfpaReduction/kfpaDataReduceGuide- 
HDec01.pdf 

3 http://www.radionet-eu.org/mwild/ParselTongue 

4 http://www.exelisvis.com/ProductsServices/IDL.aspx 



4 L. K. Morgan et al. 




Figure 1. Extracts from the 850 |xm submillimetre SCUBA map from lHatchell et all 120051) overlaid with circles showing the coverage of each observed map 
and contours showing integrated intensity of NH3 at levels of 50, 80, 90, 97 and 99% of the individual map peak. 
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Table 1. Pointing centres and observation dates for each mapped region. 





Pointing Centre Position 


Observation 


Total Integration 


Observation 


R.A.(J2000) 


Dec.(J2000) 


Date 


Time (m) 


Perseus 


Map 01 


03:25:48.8 


+30:42:24 


20 th December 2010 


21 


Man 09 


03-27-40 


+30:12:13 


7 th March 201 1 


21 


Man rn 




+31'11'09 


6 th Tanimrv 901 1 

U J ill III ill V Z.VJ 1 1 


41 


Man 04 


03-28-42 6 


+31 06-13 


9 1 st Deremher 90 1 

£d J J_^CUCJJJL7CJ i. V./ 1 V_/ 


21 


Map 05 


03:29:03.4 


+31:14:58 


20 th December 2010 


41 


Man Qf\ 


V/J . 1 \J.-J 


TJ1.1J.JJ 


91 st Deremher 9010 


21 


Man 07 


03-29-18 5 


+31:25:13 


6 th Tamiarv 901 1 


32 


Man OR 


03-32-17 5 


+30:49:49 


9 1 st Dpremher 90 1 

_ J l/LLLilll.'Ll £J\J J VI 


21 


Man OQ 




+31"19'51 


9S th Marrh 901 1 

Z.O .VlillLll Z.XJ I 1 


41 


Man 10 


03-33-1 5 1 


+31:07:04 


28 lh March 901 1 

_o ivjuiuii _v/i i 


41 


Man 1 1 

IVJtllJ 1 1 


03-41 -46 


+31:57:22 


5 th Tanuarv 201 1 

*J J dllUCU. V — * ' 1 1 


41 


Man 19 
ivi au i z. 




+39-04-06 


7 th Marrh 901 1 

/ IVlCllUil Z.VJ 1 1 


41 


Man 13 


03-44-48 8 


+32 00-29 


7 th March 201 1 

/ Aldl L 11 iUl 1 


41 


W3 










Map 01 


02:19:53.1 


+61:01:55 


16 th December 2010 


41 


Map 02 


02:20:41.0 


+61:09:42 


14 th January 201 1 


41 


Map 03 


02:21:04.0 


+61:06:01 


14 th January 201 1 


41 


Map 04 


02:21:06.0 


+61:27:28 


20 th December 2010 


41 


Map 05 


02:21:41.0 


+61:05:44 


14 th January 201 1 


41 


Map 06 


02:22:23.0 


+61:06:12 


14 th January 201 1 


41 


Map 07 


02:23:28.6 


+61:12:04 


21 st December 2010 


41 


Map 08 


02:25:31.2 


+62:06:20 


16 th December 2010 


82 


Map 09 


02:25:37.7 


+61:13:51 


19 th December 2010 


41 


Map 10 


02:26:40.0 


+62:07:00 


29 th March 2011 


41 


Map 11 


02:26:44.6 


+61:29:44 


20 th December 2010 


41 


Map 12 


02:27:17.8 


+61:57:14 


21 s ' December 2010 


41 


Map 13 


02:28:00.0 


+61:24:00 


29 th March 2011 


41 


Map 14 


02:28:12.3 


+61:29:40 


19 th December 2010 


41 


Map 15 


02:28:26.2 


+61:32:14 


20 th December 2010 


41 



3.1 Source Extraction 

In order to locate individual NH3 sources, the Starlink -CUPIljj 
implementation of the Clumpfind dWilliams et al J 1 19941) package 
was run on integrated emission maps produced from the reduced 
spectral cubes. The reasons for running the two-dimensional im- 
plementation of Clumpfind rather than a full three-dimensional 
analysis on our spectral cubes are two-fold; firstly, some results 
show that the three-dimensional implementation of Clumpfind 
may underestimate the mass of clumps in star formation regions 
dSmith et"al". 2008) and are not as robust as the two-dimensional 
version dPineda et alj|2 009). particularly in cases of limiting reso- 
lution (c.f. lWard et al.l ( 120121) . who find that the three-dimensional 
extraction is the most accurate implementation when recovering 
simulated clumps). Also, we found that no clump velocity sep- 
arations exceeded the maximum linewidth found in the relevant 
map from which the clumps were extracted. Thus, the separations 
of any multiple spectral components along the line-of-sight, are 
smaller than or approx imately equal to typical clump linewidths. 
Rosolowskv et al I l l2008l) performed a single-pointing survey of 
suspected cores in Perseus. Of the 193 sources in their catalogue, 
only six showed unambiguous evidence of multiple components 
with a mean velocity separation of 0.63 km/s. The mean observed 
linewidth of sources in our sample of Perseus cores is 0.64 km/s 
(see Tabled and so we are unlikely to be able to separate any mul- 
tiple components in this region, though they are likely present. 



5 http://docs.jach.hawaii.edu/star/sun255.htx/sun255.html 



Clumpfind was run using a lowest contour level of three times 
the r.m.s. off-source noise (cr) measurements within the completely 
sampled portion of each individual map and successive contours 
spaced by lcr. The minimum pixel count required to define a source 
was 7 (6") pixels. Where individual maps overlap, sources in the 
overlap region were extracted from the mosaicked image rather 
than the individual Daisy scans. Each extracted source was exam- 
ined in order to identify and discard any sources unduly affected by 
the higher noise values at the edges of our maps. 



3.2 Source Properties 

With the above search parameters, Clumpfind detected 84 dis- 
crete ammonia sources in Perseus and 54 in W3. Applying the 
masks produced by Clumpfind to the parameter maps which re- 
sulted from our fitting analysis routines (Section I2.2t allows us 
to extract the distribution of column density, kinetic temperature, 
etc. for each detected core. In this paper we present the integrated 
properties of these extracted cores and leave analysis of any spa- 
tial variation of these parameters within the cores for subsequent 
publications. The central positions of each source, measured in- 
tegrated intensities, fitted linewidths LSR veloc ities and any rele- 
vant m atches from the point source catalogue of iRosolowskv et al] 
d2008l) are presented in Tables [2] and [3] The Perseus and W3 re- 
gi ons have previously be en o bserved at submillim etre wavelengths 
bv lHatchelleTai] ( l2005h and iMoore et al. I J2007h respectively, the 
sources identified by those authors were cross-checked for asso- 
ciation with the ammonia sources presented in this paper. Where 
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Figure 2. The 850 |xm submillimetre SCUBA map from lMoore et alj 120071) overlaid with circles showing the coverage of each observed map and contours 
showing integrated intensity of NH3 at levels of 50, 80, 90, 97 and 99% of the individual map peak. 



a sour ce position listed in th e catalogues of either lHatchell et al] 
< l2005t) or iMoore et al.l j2007h is within the lowest contour of the 
ammonia emission as extracted by Clumpfind in this work, we 
have also listed that association 

Derived physical properties, integrated over the NH3 source 
extent are presented in Tables [4] and [5] Turbulent veloc- 
ity dispersions of each source were estimated by decon- 
volving the thermal contributions to the observed linewidths, 
dVthcrm = -^8k B r k ln2/m N H 3 , where 7\ is the mean kinetic temper- 
ature averaged over the source in question and m NH , is the mean 
molecular mass of an ammonia molecule (17.03 amu). A statisti- 
cal summary of the NH 3 source properties is given in Tableland 
their distributions are presented in Figure [4] The properties of the 
two samples appear distinct and this is confirmed by KS tests indi- 



cating that the physical parameters of the Perseus and W3 samples 
represent distinct populations with confidence levels of 5cr, except 
for column density, which is distinct with a confidence level of 4<x 

In general, the Perseus sources have smaller linewidths and 
turbulent velocity dispersions (cr tulb ), higher optical depths and fill- 
ing factors than those in W3. These results might be expected for 
any two given star-forming regions at the relative distances of W3 
and Perseus. The observations of Perseus map a physical scale 
a factor of 7.7 smaller than those of W3, assuming distances to 
Perseus and W3 of 260 pc and 2 kpc respectively, and it is to be ex- 
pected that there will be systematic resolution effects upon the mea- 
sured core properties. For instance, filling factors may be expected 
to be lower in more distant sources where there is unresolved struc- 
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Table 2. Central positions and matching sources for each mapped source in Perseus. 



Central Position 


W(NH 3 (1,1)) 


dV 


^LSR 


Matching Submm 


Matching Source 


Source R.A.(J2000) Dec.(J2000) 


(K km/s) 


(km/s) 


(km/s) 


Source (Hatchell et al. 2007) 


(Rosolowskv et al. 2008) 



P01 


03-25 29 5 


30:44:47 


7.92+0.52 


0.78 


4.11 


31 




P02 


03:25:36.9 


30:45:23 


2.95±0.52 


0.73 


4.32 


27, 28 


R10, R12 


P03 


03:25:39.8 


30:43:53 


7.82±0.52 


0.75 


4.73 


29, Bololl 


R15, R14 


P04 


03-25-48 1 


30-45-25 


2.47+0.52 


0.61 


4.45 






P05 


03:25:50.3 


30:42:32 


6 00+0 52 


0.59 


4.28 


32 


R17, R18, R19 


P06 


03:25:55.7 


30:41:02 


1.86+0.52 


0.81 


3.99 




R20 


P07 


03:26:00.3 


30:42:07 


2.12±0.52 


0.95 


4.21 






P08 


03 26 00 6 


30:40:51 


1 85+0 52 










P09 


03:27:31.0 


30:15:09 


4.42±0.51 


0.42 


4.87 




R29 


P10 


03-27-35 


30:10:22 


7 68+0 51 


0.52 


5.02 




R30 


PI 1 


03-27-36 5 


30- 13-58 


3 84+0 5 1 

J.uTXV/.J 1 


0.81 


4.55 


39 


R3 1 


P12 


03-27-37 8 


30 08-37 


3 95+0 51 


0.53 


5.28 






P13 


03:27:41.3 


30:12:17 


8.31+0.51 


0.70 


4.79 


35 36 40 


R32, R33, R34 


P14 


03-27-45 7 


30:08:39 


2.21+0.5 1 


0.73 


5.30 






P15 


03:27:49.7 


30:11:40 


6 07+0 5 1 


0.66 


4.82 


37 


R35 


P16 


03-27 54 6 


30:09:46 


2.45±0.51 


0.65 


4.60 






P17 


03:27:56.6 


30:11:15 


4.47+0.51 


0.51 


4.89 






P18 


03:27:56.8 


30:10:24 


3 gn+0 5 1 


0.60 


4.64 






P19 


03:28:28.6 


31:10:07 


j 85+0 37 


0.34 


6.98 






P20 


03:28:30.7 


31:15:14 


3.15+0.37 


0.40 


7.40 






P21 


03:28:31.8 


31:11:26 


3.34+0.37 


0.42 


7.02 


74 


R40 


P22 


03:28:33.7 


31:13:12 


3.71+0.37 


0.47 


7.13 


49 


R44 


P23 


03:28:33.8 


31:04:39 


3 67+0 37 


0.46 


6.57 


Bolo26 


R41 


P24 


03:28:35.0 


31:06:47 


1 80+0 37 


0.44 


6.67 


69 


R43 


P25 


03:28:39.5 


31:05:50 


5 40+0 37 


0.49 


6.87 


71 


R46 


P26 


03:28:44.3 


31:06:10 


4.18+0.37 


0.59 


7.1 1 


75 


R50 


P27 


03:28:49.9 


31:12:50 


2.32±0.37 


1.28 


7.16 






P28 


03:28:52.6 


31:14:51 


7.21±0.37 


1.04 


7.54 


44 


R58 


P29 


03:28:56.5 


31:13:57 


6.07±0.37 


0.83 


6.98 






P30 


03:29:01.7 


31:12:23 


4.48+0.37 


0.72 


6.95 


65 


R65 


P31 


03:29:02.6 


31:15:26 


9.58±0.37 


1.22 


7.53 


43, 52 


R67, R68 


P32 


03:29:08.3 


31:17:22 


4.95±0.37 


0.72 


8.16 


46, 62, Bolo44 


R72 


P33 


03 29 09 


31:15:21 


2 05+0 37 


1.02 


7.63 


50, 51 


R70, R73 


P34 


03:29:13.4 


31:13:20 


6 32+0 37 


0.98 


7.35 


41, 42, 48, 59, 72 


R75, R78, R80, R83 


P35 


03:29:17.1 


31:28:15 


4 49+0.36 


0.40 


7.34 


61 


R81 


P36 


03:29:17.7 


31:25:51 


3 01+0 36 


0.55 


7.1 1 


57 


R82 


P37 


03-29- 19 5 


31 -23-36 

J 1 .JU 


1 97+0.36 


0.59 


7.35 


63, 67 


R84 


P38 


03-29-20 9 


31:29:22 


2 1 0+0 36 










P39 


03-29-22 6 


31:26:11 


2.72±0.36 


0.35 


7.17 






P40 


03-29 26 


31:28:14 


3 53+0 36 


0.40 


7.30 


64 


R88 


P41 


03:32:16.1 


30:49:38 


5.52±0.47 


0.55 


6.63 


76 


R103 


P42 


03:32:55.3 


31:21:02 


3.10±0.83 


0.72 


6.34 






P43 


03:32:57.8 


31:04:10 


5 79+0 86 


0.40 


6.42 






P44 


03:32:58.6 


31:20:43 


2.90±0.83 


0.70 


6.37 




R112 


P45 


03-32 59 9 


31:03:36 


6 1 0+0 86 


0.40 


6.45 






P46 


03-33 02 2 


31:20:13 


2 84+0 83 


0.77 


6.44 






P47 


03:33:04.7 


31:05:01 


6.93±0.86 


0.50 


6.41 


5 


R113, R114 


P48 


03-33 05 7 


31:06:28 


4 32+0 86 


0.48 


6.32 




Rl 15 


P49 


03:33:10.5 


31:19:47 


5.17±0.83 


0.50 


6.62 


82 


R118 


P50 


03:33:14.5 


31:07:12 


0.28±0.86 


0.64 


6.24 


4,7 


R119 


P51 


03:33:17.4 


31:09:20 


7.72±0.86 


0.72 


6.21 


1 


R121 


P52 


03:33:21.0 


31:07:37 


3.54±0.86 


0.74 


6.45 


2 


R123 


P53 


03:33:24.3 


31:05:55 


4.58±0.86 


0.48 


6.59 




R124 


P54 


03:33:25.0 


31:19:51 


2.95±0.83 


0.89 


6.29 




R125 


P55 


03:33:26.3 


31:06:59 


3.95±0.86 


0.69 


6.61 


10 


R126 


P56 


03:33:31.0 


31:20:22 


3.58±0.83 


0.82 


6.30 




R127 


P57 


03:41:29.8 


31:57:42 


1.83±0.42 










P58 


03:41:32.9 


31:59:09 


1.61±0.42 










P59 


03:41:33.5 


31:58:05 


1.48±0.42 










P60 


03:41:36.4 


31:54:41 


1.48±0.42 










P61 


03:41:36.5 


32:00:07 


1.54±0.42 










P62 


03:41:41.1 


31:57:58 


3.64±0.42 


0.32 


9.23 


Bolo92 


R145 


P63 


03:41:46.8 


31:57:24 


3.62±0.42 


0.42 


9.24 


Bolo94 


R147 
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■l 
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Figure 3. cont. The Perseus and W3 integrated intensity maps. 



ture, while linewidths may be larger due to the inclusion of larger 
physical scales. 



3.3 Source Associations 

In the Perseus region, 6 submillimetre sources from the catalogue 
of lHatchell et alj d2007t) are covered by our observations and all of 
these are associated with an ammonia source. These sources are 
evenly split between individual matches and multiple associations, 
with 30 submillimetre sources matching a single ammonia source 



Dust and Gas in Star-Forming Environments 1 1 



Table 2. (cont.) Central positions and matching sources for each mapped source in Perseus. 





Central Position 


W(NH 3 (1,1)) 


dV 


^LSR 


Matching Submm 


Matching Source 


Source 


R.A.(J2000) 


Dec.(J2000) 


(K km/s) 


(km/s) 


(km/s) 


Source (Hatchell et al. 2007) 


(Rosolowskv et al. 2008) 


P64 


03:41:48.0 


31:54:20 


1.76±0.42 


_ 


_ 






P65 


03:41:53.2 


32:00:25 


1.89±0.42 


_ 


_ 






P66 


03:41:53.4 


31:54:33 


1.68±0.42 












P67 


03:41:55.3 


31:59:24 


1.94±0.42 


0.42 


9.20 






P68 


03:41:58.8 


31:58:54 


2.36±0.42 


0.54 


9.13 




R148 


P69 


03:42:00.4 


31:57:47 


1.73±0.42 


0.61 


9.07 






P70 


03:43:39.4 


32:03:10 


4.09±0.48 


0.80 


8.50 


23 


R156 


P71 


03:43:48.0 


32:03:18 


5.50±0.48 


0.62 


8.49 


15, 24, 26 


R157, 158, 159, 160 


P72 


03:43:55.9 


32:00:45 


4.24±0.48 


0.78 


8.69 


12 


R161 


P73 


03:43:57.3 


32:04:03 


3.76±0.48 


0.75 


8.26 


17 


R163 


P74 


03:44:01.0 


32:02:43 


5.02±0.48 


0.62 


8.47 


13, 16, 18, 21 


R162, 164, 165 


P75 


03:44:05.3 


32:01:34 


5.73±0.48 


0.62 


8.40 


20 


R169 


P76 


03:44:31.0 


32:00:35 


1.68±0.45 


0.66 


9.00 






P77 


03:44:36.9 


32:00:59 


2.00±0.45 


0.53 


9.11 






P78 


03:44:37.4 


31:58:42 


3.00±0.45 


0.68 


9.71 


19 


R176 


P79 


03:44:39.8 


31:59:36 


1.73±0.45 


0.90 


9.55 






P80 


03:44:44.1 


32:01:17 


2.03±0.45 


1.08 


9.57 


14 


R178 


P81 


03:44:49.4 


32:00:40 


2.45±0.45 


0.50 


8.79 


25 


R180 


P82 


03:44:56.1 


32:00:25 


2.35±0.45 


0.53 


8.88 




R181 


P83 


03:45:01.1 


32:00:42 


2.95±0.45 


0.40 


8.93 






P84 


03:45:06.3 


32:00:35 


2.70±0.45 


0.54 


8.93 







and 30 sources matched to 11 ammonia sources. In addition, we 
find 43 ammonia sources that d o not have an associa ted submil- 
limetre source in the catalogue of lHatchell eTai] | |2007|) . Twelve of 
these sources are weakly detected and have no (2,2) line detection 
and so the lack of submillimetre detection is unsurprising. Inspec- 
ti on of the position of th e remaining detections in the original maps 
of lHatchell elatl J2005h shows that weak submillimetre emission is 
often present below the detection limits used in the creation of the 
catalogue of lHatchell et alj d2007t) . In a few cases there is no ap- 
parent submillimetre emission whatsoever, although in these cases 
there is invariably nearby 'bowling' associated with nearby bright 
sources due to the observing techniques used in the collection of 
the data (i.e. chopping artefacts). 

The ammonia maps of the W3 region cover the positions of 
19 2 of the 316 submi llimetre sources presented in the catalogue 
of iMoore et all d2007l) . 77 of which are detected by us in am- 
monia emission, leaving 115 not detected. Compared to Perseus, 
there is a higher proportion of multiple source associations. There 
are 25 one-to-one matches and 52 submillimetre sources which 
are matched with 17 NH3 sources. Twelve sources w ere found in 
ammo nia that were not detected in the catalogue of IMoore et alj 
d2007t) , of these, 8 are weak with no (2,2) detection and 2 are close 
to the edge of the submillimetre map and so there are reasonable ex- 
planations for the non-detection of a submillimetre source. The re- 
maining two sources (W26 and W47) are nearby bright submillime- 
tre sou rces (sources 109 and 213 in the catalogue of IMoore et alj 
d2007t) respectively) which they may possibly be loosely associ- 
ated with. These submillimetre sources are more strongly associ- 
ated (closer, with more morphological similarity) to two other am- 
monia sources (W23 and W46). There is clear separation of the am- 
monia source pairs in each case and so these sources have not been 
amalgamated. There is a single case in our entire set of observations 
in which two ammonia sources appear to be clearly associated with 
a single submillimetre source. The sources W35 and W38 are both 
associ ated with the source 175 from the catalogue of IMoore et alj 
d2007h . 

In Section|2]it was stated that the coverage of sources in both 



the Perseus and W3 regions was intended to be mainly random, 
with some bias towards clustered regions. In order to test for possi- 
ble biases in our source selection we have plotted the distributions 
of submillimetre integrated fluxes for both regions in Figure [5] In 
terms of the sources covered by our observations and those that are 
not, it can be seen that a high proportion of sources may be missing 
from our observations at the low flux end. However, comparing the 
sources covered by our survey in Perseus with those in W3, we see 
no obvious biases in the initial source selection which might lead 
to, e.g. a larger average source size in W3 or similar. 

A histogram comparison between the integrated submillimetre 
fluxes of sources detected (red) and not-detected (blue) in ammo- 
nia in this region is shown in Figure [6] (no corresponding plot is 
presented for Perseus as all observed submillimetre sources within 
that region were detected). The plot shows that it is generally the 
sources weaker in submillimetre continuum which are not detected 
by this survey although there is considerable overlap between the 
ammonia detected and non-detected sources and so the separation 
of these sources may not be completely explained by a NH 3 detec- 
tion limit. 

Figure [7] shows the relation of ammonia column density to 
integrated intensity with reference to the submillimetre fluxes of 
each source, where an association has been made. The overlaid 
lines of best fit show the overall correlation of column density and 
integrated intensity in both the W3 and Perseus regions. Where 
a submillimetre association has been identified for an ammonia 
source the source is plotted as a circle, with the size of the cir- 
cle representing the 850 \vca integrated flux (for multiple associa- 
tions, the summed 850 yum integrated flux). The limiting values for 
each sample are plotted as circles with dashed lines delineating the 
lower values associated with our detecti on limits. The lower l imits 
for detection at 850 \xm are taken from lHatchell et aL and 
IMoore et al. I d2007l) and the lower limits for the detection of am- 
monia integrated intensity represent three times the median r.m.s. 
values in the integrated intensity images for the relevant region, the 
value used as a minimum threshold in the Clump find based ex- 
tractions. As there is no unique lower limit on column density mea- 
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Table 3. Central positions and matching sources for each mapped source in W3. 





Central Position 


W(NHi(l,l)) 


dV 


Vlsr 


Matching Submm 


Source 


R.A.(J2000) 


Dec.(J2000) 


(K km/s) 


(km/s) 


(km/s) 


Source (Moore et al. 2007) 


W01 


02-20-34 9 


61:26:43 


1.26±0.27 






27 


W02 


02:20:39.9 


61:09:59 


0.82±0.18 


0.67 


-48.93 


29 


W03 


02:20:44.2 


61:27:06 


1.06±0.27 


1.03 


-51.92 


30 


W04 


02:20:44.6 


61:25:58 


1.02±0.27 








W05 


02-20-48 6 


61:12:24 


0.72±0.18 








W06 


02-20-53 1 


61:26:58 


1.56±0.27 


1.75 


-51.05 


31 


W07 


02:21:00.9 


61:26:58 


1.59±0.27 


1.43 


-51.48 


35 


W08 


02:21:04.4 


61:27:41 


1.38±0.27 


0.97 


-51.17 


38 


W09 


02:21:05.8 


61:05:58 


1.12±0.18 


0.78 


-49.86 


37 


W10 


02:21:35.8 


61:05:33 


1.30+0.18 


0.96 


-49.97 


39, 42, 43, 44, 45 


Wll 


02:21:53.9 


61 -06-23 


0.96±0.18 


1.20 


-50.01 


47, 49, 50 


W12 


02:22:03.4 


61:07:18 


0.66±0.18 






51 


W13 


02:22:24.7 


61:06:04 


0.70±0.18 


0.82 


-49.99 


54 


W14 


02-23-20 6 


61:12:41 


1.13±0.19 


0.58 


-50.11 


60 


W15 


02:23:27.2 


61:12:11 


0.96±0.19 


0.67 


-50.11 


61 


W16 


02-24-58 6 


62:05:07 


1.73+0.37 


0.97 


-39.68 




W17 


02:25:07.3 


62:05:29 


1.25±0.37 


1.05 


-39.92 


78 


W18 


02:25:25.7 


61:15:05 


1.63+0.21 


1.04 


-48.86 


85, 86, 91 


W19 


02:25:27.6 


61:16:25 


1.33+0.21 






87, 92, 98 


W20 


02:25:29.2 


62:07:23 


1.30±0.37 








W21 


02-25-29 7 


61:11:01 


1 .10+0.21 


1.10 


-48.73 


90 


W22 


02:25:31.1 


61:13:03 


1.71+0.21 


1.18 


-48.30 


100, 95, 97 


W23 


02:25:32.2 


62:06:03 


3.12±0.37 


3.50 


-42.84 


109, 99 


W24 


02:25:40.4 


61:10:29 


1.50+0.21 


1.22 


-47.92 




W25 


02:25:40.7 


61:13:55 


2.02±0.21 


0.92 


-47.83 


103, 105, 107, 111, 112 


W26 


02:25:51.4 


62:06:18 


1.36±0.37 


2.08 


-38.77 




W27 


02:25:52.9 


62:04:17 


2.93±0.37 


1.75 


-38.78 


1 19 


W28 


02:26:02.2 


62:08:43 


1.97 ±0.37 


1.41 


-38.42 


127 


W29 


02:26:20.8 


62:04:21 


1.99±0.37 


2.06 


-39.84 


142 


W30 


02-26-26 3 


62:05:18 


1.77 ±0.37 


1.64 


-39.38 


148 


W31 


02:26:30.1 


61:29:38 


1.20±0.34 


0.85 


-47.08 


150 


W32 


02:26:33.3 


61:26:48 


1.42±0.34 








W33 


02:26:35.6 


62:06:42 


1.51 ±0.37 


1.42 


-39.15 


152, 161, 166 


W34 


02:26:37.3 


61:25:46 


1.55±0.34 








W35 


02:26:41.6 


61:32:49 


1.31+0.34 


0.94 


-47.88 


175 


W36 


02:26:42.1 


61:25:23 


1.77±0.34 








W37 


02:26:43.1 


61:33:34 


1.23±0.34 






171 


W38 


02:26:43.4 


61:31:49 


1.17±0.34 


1.11 


-47.79 


175 


W39 


02:26:44.2 


61:29:48 


1.80±0.34 


1.25 


-47.69 


172 


W40 


02:26:45.8 


62:08:41 


1.26±0.37 


1.41 


-38.68 


176, 180 


W41 




61 '25'39 


1.62±0.34 








W42 


02:26:49.8 


62:10:04 


1 40+0.37 








W43 


02:26:59.3 


61:54:45 


1.93±0.47 


0.80 


-45.45 


196 


W44 


02:26:59.6 


61:29:51 


2.03±0.34 


1.78 


-47.69 


195,214 


W45 


02:26:60.0 


61:53:50 


1.68±0.47 


1.99 


-45.17 


198 


W46 


02:27:02.6 


61:52:03 


4.25±0.47 


2.78 


-47.20 


213 


W47 


02:27:11.2 


61:53:05 


3.24±0.47 


2.47 


-48.48 




W48 


02:28:01.8 


61:24:34 


2.29±0.37 


1.84 


-48.05 


262, 270, 275,280, 281,286 


W49 


02:28:01.9 


61:27:22 


5.22±0.37 


1.61 


-48.63 


266, 268, 284, 288 


W50 


02:28:05.3 


61:29:58 


4.49±0.37 


1.86 


-47.85 


264, 285, 287, 290, 292 


W51 


02:28:09.8 


61:25:41 


2.84±0.37 


2.11 


-47.49 


289, 293 


W52 


02:28:15.1 


61:26:38 


2.83±0.37 


2.46 


-48.02 


291 


W53 


02:28:23.1 


61:31:52 


2.19±0.37 


1.94 


-47.04 


297, 299, 300 


W54 


02:28:58.7 


61:33:29 


1.11±0.19 


0.81 


-51.46 


315 



surements, the plotted limit is simply taken as the lowest measured 
value in the relevant region. 

A correlation is found between the 850 |J.m integrated flux and 
the ammonia integrated intensity for both the W3 and Perseus re- 
gions, with slightly differing slopes in these correlations for W3 
and Perseus, due to the higher average temperature of cores in W3. 
Where no associated submillimetre source was found, the ammo- 



nia source is plotted as a star (W3) or diamond (Perseus). The plot 
shows that, generally, the sources which are detected in ammonia 
but do not have an 850 (J.m assocation are lower in integrated inten- 
sity and column density. There is considerable overlap between the 
subsamples of submillimetre associations/non-assocations but the 
overall trend of submillimetre assocations representing the highest 
intensity/column density sources is clear. 
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Table 4. Derived physical properties for each mapped source in Perseus. 





7k 


Optical 


Filling 


Column Density 


Source 


(K) 


Depth 


Factor 


(/10 13 cm" 2 ) 


P01 


11.30 


4.13 


0.70 


42.87 


P02 


11.37 


4.01 


0.81 


46.44 


P03 


11.40 


2.85 


0.60 


36.05 


P04 


13.47 


1.64 


0.22 


22.02 


P05 


10.19 


4.14 


0.62 


34.65 


P06 


15.69 


1.21 


0.21 


18.16 


P07 


18.56 


1.75 


0.07 


54.21 


P08 


_ 




. 


m_ 


P09 


11.39 


3.04 


0.47 


20.81 


P10 


12.88 


1.16 


0.49 


10.69 


PI 1 


13.36 


1.60 


0.38 


25.35 


P12 


12.92 


1.81 


0.53 


15.87 


PI 3 


12.69 


3.10 


0.60 


35.74 


P14 


19.99 


0.93 


0.26 


21.65 


P15 


12.75 


2.19 


0.61 


22.89 


P16 


17.26 


0.87 


0.33 


18.94 


PI 7 


11.84 


2.29 


0.64 


18.33 


P18 


13.74 


1.54 


0.56 


15.92 


P19 


10.35 


3.66 


0.35 


17.79 


P20 


10.84 


5.49 


0.37 


34.49 


P2I 


11.58 


3.50 


0.32 


23.10 


P22 


12.00 


3.37 


0.34 


27.37 


P23 


11.41 


3.34 


0.42 


25.25 


P24 


11.61 


4.40 


0.15 


30.65 


P25 


10.45 


3.62 


0.53 


27.23 


P26 


11.82 


3.76 


0.41 


32.52 


P27 


17.82 


1.10 


0.29 


43.15 


P28 


14.69 


1.57 


0.52 


38.32 


P29 


14.17 


1.75 


0.48 


28.66 


P30 


13.89 


1.85 


0.44 


23.41 


P31 


16.45 


1.93 


0.45 


55.50 


P32 


14.56 


2.08 


0.46 


30.47 


P33 


13.74 


1.95 


0.75 


42.11 


P34 


13.19 


1.49 


0.59 


30.53 


P35 


12.03 


3.21 


0.46 


22.51 


P36 


12.63 


2.10 


0.36 


21.50 


P37 


14.86 


1.97 


0.14 


24.40 


P38 


^ _ 


_ 


_ 


_ 


P39 


12.63 


1.98 


0.34 


12.98 


P40 


10.99 


3.11 


0.45 


21.72 


P4I 


12.07 


2.86 


0.51 


31.34 


P42 


17.44 


0.48 


0.45 


11.83 


P43 


10.40 


4.54 


0.54 


27.49 


P44 


14.09 


0.82 


0.51 


11.06 


P45 


10.01 


4.36 


0.61 


27.05 


P46 


12.55 


1.48 


0.28 


20.63 


P47 


10.30 


4.57 


0.60 


32.43 


P48 


10.82 


2.66 


0.48 


19.73 


P49 


10.14 


5.15 


0.42 


33.99 


P50 


11.63 


4.50 


0.62 


44.07 


P51 


12.22 


2.70 


0.53 


34.22 


P52 


11.18 


4.44 


0.72 


52.61 


P53 


10.52 


3.99 


0.45 


26.56 


P54 


12.38 


1.76 


0.25 


28.44 


P55 


12.41 


1.54 


0.48 


19.14 


P56 


12.88 


1.96 


0.26 


33.21 


P57 










P58 










P59 










P60 










P6I 










P62 


9.45 


4.78 


0.42 


22.59 


P63 


10.49 


3.35 


0.39 


22.87 
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Table 4. (cont.) Derived physical properties for each mapped source in Perseus. 





7k 


Optical 


Filling 


Column Density 


Source 


(K) 


Depth 


Factor 


(/10 13 cm" 2 ) 


P64 




_ 




_ 


P65 


_ 


_ 







P66 














P67 


14.33 


1.97 


0.11 


15.37 


P68 


14.52 


1.00 


0.25 


8.50 


P69 


15.39 


1.31 


0.24 


9.59 


P70 


16.52 


1.32 


0.36 


23.58 


P7I 


13.49 


2.28 


0.48 


22.87 


P72 


16.38 


1.99 


0.32 


34.25 


P73 


15.18 


1.71 


0.29 


25.46 


P74 


13.25 


2.27 


0.44 


23.89 


P75 


12.16 


2.66 


0.47 


31.10 


P76 


16.88 


0.78 


0.21 


10.81 


P77 


13.59 


1.02 


0.30 


10.76 


P78 


13.23 


1.55 


0.35 


19.04 


P79 


18.86 


0.53 


0.20 


15.41 


P80 


15.37 


0.61 


0.24 


17.85 


P81 


12.62 


1.90 


0.28 


17.64 


P82 


12.13 


1.21 


0.43 


9.32 


P83 


11.76 


2.17 


0.42 


14.60 


P84 


13.21 


2.66 


0.22 


24.50 



Comparing the submillimetre fluxes of the sources found in 
W3 which do not have detectable ammonia association (no such 
sources exist in our Perseus samples) shows that it is the weaker, 
lower mass sources which are not detected in ammonia. The in- 
tegrated intensity completeness limits of our respective ammonia 
samples are estimated from the turnover in the integrated inten- 
sity number distributions at 1.6 K km/s for W3 and 2.5 K km/s for 
Perseus, corresponding to logi integrated intensity values of 0.2 
and 0.4 respectively on the plot presented in Figure [7] The differ- 
ing completeness limits reflect the higher noise found in our ob- 
servations of the Perseus region. From the plot we can see that 
few of our Perseus sources lay below the Perseus completeness 
limit while there is a significant proportion of W3 sources which 
fall the W3 completeness limit. This difference in the proportion 
of sources above the respective completenes s limits expla i ns the 
number of sources detected at 850 \im by iMoore et alj ( 2007) 
which are not detected in this study. With the correlations found 
between the integrated ammonia intensity and 850 \im integrated 
flux measurements the completeness of our overall sample is lim- 
ited by a combination of the relevant completeness limits of our 
ammonia observations and the relevant submillimetre observations. 
From a comparison of the estimated completeness limit in our am- 
monia source sam ple t o the complet e ness limits established by 
IMoore et al. | j2007l) and lHatchell et all J2007I) we have determined 
that the completeness of our sample is limited by submillimetre 
sensitivity in Perseus and ammonia sensitivity in W3. 

3.4 The Relationship of Submillimetre and Ammonia 
Emission 

The offsets between the position of peak intensity of each ammonia 
source and that of its associated submillimetre source are presented 
in Figure [8] in the case of multiple submillimetre source associa- 
tions, the closest source is used. The plot shows that, while sources 
share common regions of emission between the submillimetre and 
ammonia, peaks of emission may be separated by some distance. 
Approximately 60% of the sources are within a radius equivalent 



to that of the GBT beam. The distributions of peak offsets in the 
Perseus and W3 regions are plotted in Figure [9] It can be seen 
that peak offsets are skewed towards lower values overall, though 
peak offsets clearly do not peak at zero in either region. Thus, the 
peak positions of submillimetre emission do not appear tightly con- 
strained to the position of ammonia emission. 

The relative distributions of submillimetre and NH 3 emission 
over six sources, three from Perseus (Perseus 05, 33 & 41) and 
three from W3 (W3 18, 25 & 49), are shown in Figure[l0]with sub- 
millimetre emission contours overlaid upon NH 3 intensity images. 
These sources are representative of each region in that the ammonia 
emission associated with Perseus sources is typically significantly 
larger than the extent of submillimetre emission (taking the 50% 
surface brightness contour as our size metric), for W3 sources there 
is typically less of a disparity. The sources Perseus 33 and all of the 
W3 images in Figure[l0]show the presence of multiple submillime- 
tre cores encompassed within the 50% ammonia emission contour. 
Equivalent images for all associated sources are presented in Figure 
lAlin online material. 

Inspection of Figures \10\ and [A] suggests that the ammonia 
emission in Perseus cores is more spatially extended than the re- 
lated submillimetre emission while in W3 the contrast is less dis- 
tinct. Any difference between the distribution of gas and dust, or 
their tracers, is of interest to star formation studies as such trac- 
ers are used regularly, and often interchangeably, to determine the 
physical properties of protostellar cores. Therefore we wish to es- 
tablish metrics for the sizes of sources as traced by ammonia and 
submillimetre continuum and their relationship to each other. 

It has been noted by previous authors (seelReid et alj2010] and 
references therein, also lCurtis & Richer! 20 id) that different source 
extraction techniques may produce significantly different results. It 
may therefore be important to state that all sources identified here, 
in both ammonia and submillimetre emission, have been identified 
using Clump find. The size of each s ource was defined by its 50% 
surface-brightness contour. Following lHatchell et al. 1 120071) we de- 
fine R e (f, the effective radius of a source, as -./-, where A is the 
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Table 5. Derived physical properties for each mapped source in W3. 





T k 


Optical 


Filling 


Column Density 


Source 


(K) 


Depth 


Factor 


(/10 13 cm" 2 ) 


W01 


_ 


_ 


_ 


_ 


W02 


15.08 


2.70 


0.04 


41.00 


W03 


18.20 


1.65 


0.03 


71.82 


W04 












W05 














W06 


17.71 


1.56 


0.06 


84.70 


W07 


15.83 


1.02 


0.12 


33.69 


W08 


15.41 


1.29 


0.10 


32.16 


W09 


13.49 


1.42 


0.11 


23.40 


W10 


14.64 


1.29 


0.14 


25.19 


Wll 


16.12 


1.11 


0.08 


31.36 


W12 














W13 


17.40 


0.60 


0.13 


14.07 


W14 


12.44 


3.48 


0.08 


39.87 


W15 


15.10 


3.23 


0.05 


37.49 


W16 


16.91 


1.71 


0.11 


43.33 


W17 


17.24 


1.18 


0.09 


31.82 


W18 


13.75 


1.32 


0.18 


27.68 


W19 














W20 














W21 


17.60 


0.67 


0.13 


19.53 


W22 


13.88 


1.61 


0.14 


35.70 


W23 


35.04 


0.73 


0.08 


199.64 


W24 


15.66 


0.70 


0.20 


19.31 


W25 


12.54 


1.98 


0.17 


35.18 


W26 


34.73 


0.32 


0.08 


48.97 


W27 


21.67 


1.10 


0.11 


71.82 


W28 


21.88 


0.92 


0.11 


46.84 


W29 


22.83 


0.85 


0.10 


67.41 


W30 


20.74 


0.79 


0.10 


45.06 


W31 


15.10 


2.59 


0.05 


47.99 


W32 







M 





W33 


23.11 


0.49 


0.11 


27.09 


W34 













W35 


16.81 


0.94 


0.19 


21.75 


W36 













W37 














W38 


18.19 


0.43 


0.16 


12.29 


W39 


15.11 


1.05 


0.14 


28.84 


W40 


20.90 


0.40 


0.10 


20.04 


W41 













W42 


— 


— 


— 


— 


W43 


17.22 


2.26 


0.10 


43.35 


W44 


15.76 


1.10 


0.14 


40.11 


W45 


18.90 


0.30 


0.21 


18.57 


W46 


25.84 


1.39 


0.09 


189.02 


W47 


22.12 


1.37 


0.09 


134.67 


W48 


17.47 


0.94 


0.14 


46.81 


W49 


16.57 


1.27 


0.23 


50.05 


W50 


17.97 


1.13 


0.19 


56.89 


W51 


17.43 


1.05 


0.14 


64.44 


W52 


19.05 


1.13 


0.10 


84.35 


W53 


19.67 


0.78 


0.11 


54.48 


W54 


23.76 


1.18 


0.10 


35.21 



area within the 50% contour, having deconvolved the FWHM beam 
size assuming a Gaussian source profile and a Gaussian beam. This 
deconvolution method is only intended to be approximate as the 
assumption of a Gaussian source is somewhat unrealistic. 

Sizes for all ammonia sources are presented in Tables|7]and[8] 
along with the size s of the related sources f rom the catalogues of 
iMoore et alj d2007t) and lHatchell et alj d2005h . The sizes of the sub- 



millimetre sources have been re-measured from the original data 
using the same method to ensure consistenc y, also, no s i zes fo r 
the submillimetre sources were presented by IMoore et ail J2007l) . 
Where the emission of an ammonia source covers more than one 
submillimetre source, the combined area of the latter is used. Am- 
monia sizes vary surprisingly little within and between each region, 
with mean deconvolved R c( r measurements of 35.6" and 32.6" for 
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Table 6. Summary of physical properties for sources in Perseus and W3. 





Mean 


Perseus 
Median 


Std. Dev. 


Mean 


W3 
Median 


Std. Dev. 


dV (km/s) 


0.64 


0.61 


0.20 


1.44 


1.25 


0.63 


(Tturb (km/s) 


0.25 


0.24 


0.10 


0.59 


0.54 


0.27 


?k (K) 


13.14 


12.69 


2.31 


18.50 


17.43 


4.85 


Optical Depth 


2.42 


2.08 


1.24 


1.27 


1.13 


0.72 


Filling Factor 


0.42 


0.43 


0.16 


0.12 


0.11 


0.05 


Column Density (/10 13 cm -2 ) 


25.78 


23.89 


10.69 


50.07 


40.11 


39.92 




Figure 4. Histograms of the physical properties of the Perseus (red) and W3 (blue) regions. 
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Figure 5. Distributions o f the submillimetre in tegrated fluxes of the Perseus (left) and W3 (right). The solid black line shows the distribution of fluxes from the 
entire catalogue of either Hatchell e t al.l j2007t) for Perseus or iMoore et all fe007t) for W3. The solid red line shows the distribution of fluxes of those sources 
which are covered by our observations and the dashed blue line shows the difference between the two distributions. 




Figure 6. Histograms of the (log) integrated fluxes of the sources detected (red) and not-detected (blue) in the W3 region by this survey. 



the Perseus and W3 samples respectively. Median sizes agree with 
the means to significantly better than the sample standard devia- 
tion in each case, indicating that each sample is truly self-similar 
( Rrff-Rgf _ q 33 and o 05 for p erseus an(1 W3 respectively). 

Histograms of the source-size distributions, comparing the 
two tracers and the two cloud samples, are shown in Figure QT] 
The narrow range of NH 3 source sizes noted above is clearly seen. 
Also, the submillimetre and NH 3 sizes in W3 are distributed rather 
similarly whereas in Perseus, submillimetre radii are significantly 
smaller, with barely any overlap between the histograms of the 
two tracers. Kolmogorov-Smirnov (KS) tests of the submillime- 
tre/ammonia R e ff ratios suggest that the populations of W3 and 
Perseus are significantly different in terms of relative source size, 
with a confidence equivalent to 100% to six significant figures 

In order to investigate potential differences in the structure 
of the two star-forming regions on unresolved scales we com- 
pared the ratio of peak-to-integrated submillimetre fluxes for each 
source (hereafter, r] p J l " k ). A source which displays a sharply peaked 
flux profile and/or emission which is not extended relative to the 
beam should present a rf.° approaching unity. Conversely, sources 
which are large in comparison to the beam with 'flat' flux profiles 



should have lower ratios. Figure [12] shows how the values of this 
ratio are distributed for Perseus and W3 sources. It can be seen that 
sources in W3 are typically lower in rf™ k and that values in Perseus 
are both higher on average as well as more widely distributed. The 
mean value of rf^ k in Perseus is 0.28 while in W3 the value is 0.09. 
The relative distributions of rf™ k seen in this plot are very similar 
to the distributions of filling factor which can be seen in Figure [4] 
the mean values of filling factor in our sample are 0.42 in Perseus 
and 0.12 in W3. At first, low/high values of filling factor and r} p ™ k 
in W3/Perseus appear contradictory as extended sources may be 
expected to have high filling factors and low rf™ k . However, the 
correlation of filling factor and flux ratio suggest that sources in 
W3 are extended but 'lumpy', i.e. that there are multiple sources 
within a beam. This explanation is consistent with the observation 
that each ammonia source in Perseus is typically associated with 
fewer submillimetre sources than in W3. A more detailed analysis 
of source structure and variation of parameters within sources will 
be performed in a later work. 

With nearly an order of magnitude difference in the distance 
of the Perseus and W3 regions, it is possible that the contrast in the 
submillimetre source sizes and submillimetre/ammonia source size 
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Figure 7. The distribution of ammonia column density against integrated intensity on a log-log scale. Circles represent submillimetre associations, the size of 
the circle representing the log value of the 850 u-m integrated flux on a normalised scale. Blue and red circles represent sources in W3 and Perseus respectively, 
stars (W3) and diamonds (Perseus) show the sources in those regions with no submillimetre association. Lines of best fit to all data points for each region are 
overlaid. 
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Figure 8. Offsets of the peak positions of each submillimetre source in Perseus (red diamonds) and W3 (blue stars) relative to the peak ammonia position 
centred at 0,0. The circle represents the radius of the GBT beam. Typical pointing errors in our observations are 4" . 
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Figure 9. Histogram of the peak positions of submillimetre sources related to the ammonia sources in Perseus (red) and W3 (blue) 



ratios is due to the difference in resolved physical scales. In order 
to investigate this, we convolved and regridded the Perseus maps 
(both ammonia and submillimetre) to an equivalent beam and pixel 
size to place it effectively at the distance of W3 (2 kpc) and then 
re-ran the source extraction routines and subsequent analyses. The 
resulting blending of sources means that our Perseus sample size 
was drastically reduced, in almost all cases to a single source per 
KFPA map. However, the results serve as a control for the effects 
that distance may have on our study and allow us to use our Perseus 
observations as a template against which to compare our W3 obser- 
vations. 

The source sizes found for the smoothed, regridded Perseus 
maps are shown in Table [9] It can be seen that the effect of artifi- 
cially placing the Perseus sources at the distance of W3 decreases 
the measured size of sources in both submillimetre and ammonia 
emission by a factor of ~2. The submillimetre/ammonia size ratio 
increases only slightly, by less than half the standard deviation of 
the sample. This suggests that the relative distance of the Perseus 
and W3 regions is unlikely to be the cause of the difference between 
the size ratios in the two samples. 




3.5 The Fractional Abundance of Ammonia 

The fractional abundance of ammonia in protostellar clouds is a 
poorly defined quantit y with estimates in speci fic regions vary- 
ing from 7 x 10 " 10 fDi Francesco et alj l2002h to ~1 x lfr 8 
dTafalla et al]|2004 iFriesen et al . 2009). More recent values aver- 
aged over wide-field samples and Galactic plane surveys favour 
highe r abundances with averages of ~4 x 10~ 8 jPunham et alj 
l201ll) and 1.2 x 10~ 7 l lWienen et alj2012l) found in single-pointing 
ammonia observations of dust peaks found in 1 . 1 mm (the Bolocam 
Galactic Plane Survey. lAguirre et alj|201 lh and 870 \im surveys 
(the Apex Telescope Larg e Area Survey of the Galaxy, ATLAS- 
GAL, ISchuller et al1l2009h , respectively. The abundance has also 
been shown to be d ependent upon the internal structure of cores in 
some case studies dTafalla et alj|2004l) . The variation of ammonia 



abundance means that the molecule may not be a good estimator of 
total core mass. 

We have calculated the fractional abundance of ammonia for 
each source, both as an integrated property as well as within each 
source individually. The imperfect correlation of ammonia and sub- 
millimetre emission means that fractional abundances may only be 
determined for certain portions of many cores. 

The fractional abundance was calculated through the simple 
ratio of NH 3 column density to H 2 column density measured from 
submillimetre maps. The latter is determined through the equation 

N(H 2 ) = S „/ [a ml im H K v B v {T d )\ (1) 

where 5 „ represents the 850 |i.m flux density (evaluated per pixel), 
C2 m is the main beam solid angle, ji is the mean molecular weight 
(taken as 2.29), m H is the mass of a hydrogen atom, k v is the dust 
opacity per unit mass at 850 [im, taken here and elsewhere as 0.01 
cm 2 g _1 . ByiTj) represents the Planck function evaluated at a dust 
temperature T d . T d is taken to be equal to the integrated NH 3 ki- 
netic temperature over the source, this approximation is valid in the 
case of star-forming cores as, firstly, there is typically little varia- 
tion of r k over a given source and secondly, the kinetic temperature 
of ammonia in star-forming sources has be en found to be roug hly 
equivalent to the dust temperature (see, e.g. lMorgan et alfeoioh . 

It should be noted that we have used the temperature- 
averaging method described in iMorgan et all d2012l) to extend our 
maps of ammonia column density for optimum coverage of submil- 
limetre emission. We find that fractional abundances of ammonia, 
averaged over each source, range from 7.9 x 10~ 8 to 2.7 x 10~ 7 in 
our Perseus sample (mean = 1.6 x 10~ 7 ) and from 1.1 x 10~ 7 to 
3.6 x 10~ 6 in W3 (mean = 6.2 x 10~ 7 ). The distributions of abun- 
dance values for both Perseus and W3 cores are shown in Figure 
[T3l The standard deviation of abundance values in Perseus is typ- 
ically 29% of the mean while, in W3, it is 98%. In addition, the 
variation within individual cores is seen to be significantly higher 
in W3 than in Perseus with average standard deviations of 77% and 
39% of the individual mean value respectively. 

It is common practice to estimate the fractional abundance 
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Figure 12. The distributions of peak to integrated flux ratios for submillimetre sources associated with our ammonia sources in Perseus (red) and W3(blue). 




Figure 13. Histogram of the (log) fractional abundances found in Perseus (red) and W3 (blue) cores. 



of ammonia in star-forming regions through single-pointing ob- 
servations directed at previously identifi e d peaks of submil- 
limetre emission (e^g. IWienen etZI [20121 ; iDunham et alj 1201 ll ; 
iRosolowskv et alj|2008h . In order to test the assumption that mea- 
surements of ammonia emission made at the position of dust emis- 
sion peaks are representative of the average characteristics of a 



given source, we compare the abundance found at the submillime- 
tre peak position with the average value over the entire source. We 
define the quantity 5 X = [x sm - x] lx, where Xsm is the fractional 
abundance of ammonia measured at the position of peak submil- 
limetre emission and x is the fractional abundance averaged over 
the source. Figure [T4l shows the variation of 6 X . This histogram 
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Table 7. Source sizes in Perseus. 





NH 3 FU 


Submillimetre R<,fy 


Submillimetre R c (f 


Source 


(") 


(") 


NH 3 R eff 


Perseus 


POl 


27.4 


20.8 


0.76 


DA 1 

POj 


42.1 


9.5 


0.23 


P05 


40.7 


26.4 


0.65 


PI 1 


29.4 


13.5 


0.46 


P13 


45.4 


23.5 


0.52 


P15 


47.1 


20.3 


0.43 


P22 


31.8 


11.0 


0.35 


P25 


41.5 


17.4 


0.42 


P26 


31.5 


12.6 


0.40 


P2S 


29.2 


8.7 


0.30 


P30 


32.0 


9.7 


0.30 


P3 1 


41.3 


18.9 


0.46 


P32 


40.8 


20.2 


0.50 


P33 


42.4 


21.1 


0.50 


P34 


47.0 


29.4 


0.63 


P35 


41.9 


12.9 


0.31 


P36 


38.1 


16.9 


0.44 


P37 


23.5 


26.0 


1.11 


P40 


27.1 


12.8 


0.47 


P41 


27.6 


8.0 


0.29 


P47 


38.8 


37.1 


0.96 


P49 


31.6 


14.3 


0.45 


P50 


42.1 


23.4 


0.56 


P51 


30.4 


10.0 


0.33 


P52 


37.0 


13.3 


0.36 


P70 


31.0 


10.2 


0.33 


P71 


50.7 


24.0 


0.47 


P72 


29.6 


11.4 


0.39 


P73 


32.7 


16.9 


0.52 


P74 


43.0 


21.0 


0.49 


P75 


33.1 


22.9 


0.69 


P78 


24.5 


13.0 


0.53 


P80 


23.0 


9.9 


0.43 


Mean 


35.6 


17.2 


0.49 


Median 


33.1 


16.9 


0.46 


St.Dev. 


7.5 


6.9 


0.19 



shows that the value of ammonia abundance found at the position 
of peak submillimetre emission is significantly lower than the mean 
in all but one case and is not representative of the region as a whole. 
On average, the Perseus and W3 samples show similar values of S x 
with mean differences of -52% and -58% respectively. However, 
Figure [14] shows that the distributions of S x are skewed with the 
median in Perseus (6 X = -48%) higher than the mean. Conversely, 
the median average in W3 (d Y = -65% ) is low er than the mean. 
These results echo those of lFriesen et alj | |2009|) in which they find 
that dust clumps in the region of Ophiuchus, as identified by 850 
[im submillimetre emission, are associated with minima in the frac- 
tional abundance of ammonia. The possibility that this effect might 
be caused by peaks in the submillimetre smaller than the GBT beam 
was checked for and discounted by smoothing the submillimetre 
data to the resolution of the ammonia maps. 

It is particularly noteworthy that the median values of 6 X are 
comparable to the standard deviations found for fractional abun- 
dance within individual cores, indicating that the practice of taking 
the position of peak submillimetre emission as representative of an 
entire source when measuring fractional abundance will lead to a 
bias greater than 1 a in a large proportion of cases. 



3.6 Virial Ratios 

In order to further examine the nature of the dense cores in W3 in 
comparison to those in Perseus, we calculated the virial ratios of 
each source for which the necessary parameters had been obtained. 
Assuming spherical geometry we use the formulation 



to determine the virial mass of each core, where G is the gravita- 
tional constant and R c is the ammonia core radius as presented in 
Tables [7] and [8] &y m is the 3D velocity dispersion integrated over 
the core in question, calculated by substituting the thermal contri- 
bution to the measured line width from ammonia with that assum- 
ing a mean molecular weight of 2.29. Having determined the virial 
mass of each core, we then calculated the virial ratio using the ob- 
ser ved mass, de t ermin ed us ing the submill i metre fluxes published 
by iMoore et alj J2007I) and lHatchell et al] fe007l) and the kinetic 
temperatures determined from our ammonia observations, the re- 
sulting virial ratio values are plotted as a histogram in Figure [T31 
and as a function of column density in Figure [T6l 

We have already pointed out that several measures of the mass 
in a star-forming region may be found to vary between submil- 
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Table 8. Source sizes in W3. 





NH 3 R cff 


Submillimetre R<,fj 


Submillimetre 


Source 


(") 


(") 


NH 3 Reft 


W3 


W01 


26.3 


10.5 


0.40 


W03 


29.6 


19.6 


0.66 


W06 


32.7 


22 7 


0.69 


W07 


31.5 


Xl.l 


0.56 


W08 


35.7 


26.2 


0.73 


W09 


24.7 


17.4 


0.70 


W10 


27.8 


46.4 


1.67 


Wll 


35.7 


31.4 


0.88 


W12 


22.3 


20.9 


0.94 


W13 


25.2 


13.8 


0.55 


W14 


19.2 


19.4 


1.01 


W15 


21.7 


13.9 


0.64 


W17 


18.3 


46.6 


2.55 


W18 


37.6 


35.1 


0.93 


W19 


34.6 


44.0 


1.27 


W21 


41.0 


21.0 


0.51 


W22 


36.9 


38.8 


1.05 


W23 


27.6 


29.2 


1.06 


W25 


45.8 


52.2 


1.14 


W27 


22.8 


12.6 


0.55 


W28 


23.8 


26.5 


1.1 1 


W29 


25.4 


31.0 


1.22 


W30 


37.4 


26.0 


0.70 


W31 


27.4 


13.3 


0.49 


W33 


52.6 


41.8 


0.79 


W35/W38 


30.6 


52.1 


1.70 


W37 


16.4 


18.4 


1.12 


W39 


36.6 


13.6 


0.37 


W40 


36.6 


40.2 


1.10 


W43 


9.8 


16.6 


1.69 


W44 


42.2 


30.0 


0.71 


W45 


21.7 


33.1 


1.53 


W46 


36.9 


11.4 


0.31 


W48 


54.7 


53.2 


0.97 


W49 


51.7 


55.2 


1.07 


W50 


54.8 


58.2 


1.06 


W51 


32.0 


44.4 


1.39 


W52 


37.9 


38.4 


1.01 


W53 


53.5 


7.6 


0.14 


W54 


23.5 


12.6 


0.54 


Mean 


32.6 


29.1 


0.94 


Median 


32.0 


26.5 


0.94 


St.Dev. 


11.0 


14.6 


0.46 



limetre continuum and NH3 molecular line observations. For ex- 
ample, the peak positions of these two tracers are not tightly corre- 
lated in space and there is an apparent anti-correlation between the 
column density of submillimetre continuum measurements and the 
fractional abundance of ammonia. Therefore, there is some doubt 
that ammonia and submillimetre continuum emission are tracing 
the same mass component. If they are not, then use of the ob- 
served submillimetre mass to establish the virial ratio may intro- 
duce some systematic bias. We have persisted in using the submil- 
limetre mass primarily because of the lack of any more reliable 
alternative. The fractional abundance of ammonia, as found in the 
literature, is highly variable/uncertain (see section |4~3l l with values 
ranging over several orders of magnitude. Any estimate of mass de- 
rived using an assumed value of abundance is therefore subject to 
the same uncertainties. While we have far more tightly constrained 
values of abundance derived for each core in the previous section, 



these values are dependent upon the assumption that ammonia and 
submillimetre continuum emission both trace the same mass com- 
ponent and so we would have no better an estimate of the mass 
than simply using the submillimetre-derived mass. While we must 
be aware that uncertainties may be introduced to our virial ratio de- 
termination through inclusion of the submillimetre mass, analysing 
the values of M v alone show that the basic virial properties of cores 
in the Perseus and W3 regions are consistent. That is, virial masses 
in Perseus are higher than in W3. 

The results in Figures [731 and [T6l show significant and simi- 
lar scatter in the virial ratios of cores in both the Perseus and W3 
regions. This scatter may be due in part to the uncertainties in de- 
termining virial ratio values. Estimates for submillimetre-derived 
mass require estimates of k v , the frequency-dependent dust opacity 
per unit mas s. This is highly unc ertain with estim ates varying from 
0.01 cm 2 g~'( lMo"ore et alj[2007h to 0.02 cmV dThompson et alj 
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Figure 14. Histogram plot of the difference in the fractional abundance at the submillimetre peak position compared to the mean found in Perseus (red) and 
W3 (blue) cores. 
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Table 9. Convolved and regridded source sizes in Perseus 





NH 3 Reff 


Submillimetre R.ff 


Submillimetre R c ff 


Source 


(") 


(") 


NH 3 R cff 


Perseus 


MapOl 


11.9 


7.7 


0.65 


Map02 


22.8 


10.3 


0.45 


Map()7 


18.7 


7.8 


0.42 


Map08 


23.1 


6.0 


0.26 


Map09 


18.0 


16.1 


0.89 


Map 10 


10.4 


8.9 


0.86 


Map 11 


14.9 


8.5 


0.57 


Map 12 


18.3 


10.5 


0.57 


Map 13 


14.5 


8.6 


0.59 


P.b.l 


16.0 


9.6 


0.60 


P.b.2 


20.7 


8.8 


0.43 


Mean 


17.2 


9.4 


0.57 


Median 


18.0 


8.8 


0.57 


St.Dev 


4.1 


2.6 


0.19 
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Figure 16. Plot of the virial ratio vs. column density for each of our ammonia cores found in the Perseus (red diamonds) and W3 (blue stars) regions. Green 
triangles show the calculated virial ratios of our convolved Perseus sources, artificially observed at the distance of W3. The dotted line shows the line of virial 
equilibrium. 



even in similar regions. In addition to this, a significant un- 
certainty in mass can be attributed to temperature estimates through 
the Planck function. Therefore, we estimate uncertainties in virial 
ratio to be as large as 50%. However, even this uncertainty cannot 
account for the scatter seen in virial ratios in Figure [16] In both 
samples the distribution of virial ratio values is skewed similarly 
with a longer tail to lower values. However, the W3 distribution is 
shifted to higher values and a significant fraction of the cores in 
Perseus are unbound, while W3 cores are typically more strongly 
bound. Median virial ratios for each region are 0.23 and 0.72 for 
Perseus and W3 respectively, mean values are significantly higher 



at 0.53 and 1.19, reflecting the skew. In Figure [T6l we see a weak 
correlation between the virial ratio and the NH 3 column density 
in Perseus with a Pearson correlation coefficient of r=0.33 which, 
with 33 data points, is significant at around the 2-sigma level. The 
latter suggests that the higher column-density cores in Perseus are 
more dynamically similar to those in W3 than the general popula- 
tion. There is no corresponding significant correlation in the W3 
data (r=0. 13 which is not significant at the 1-sigma level). There 
is also a marginally significant anti-correlation of virial ratio with 
7\ in the Perseus sample (r=-0.29), probably reflecting the fact that 
the gas kinetic temperature is a significant fraction of the observed 



line width in those sources. No strong correlation of the virial ratio 
with any other parameter was found. 

The virial analysis of Perseus sources by I Kirk et all l l2007t) . 
using single-point N 2 H + and submillimetre continuum map data, 
found results consistent with ours, i.e. median virial ratios approxi- 
mately a factor of two less than unity. However, they find that cores 
in Perseus are typically bound when modelle d external pressure is 
taken into consideration. iRoberts et all d 19971) find virial masses in 
W3 are typically higher than measured H 2 masses by ~30%, also 
consistent with the virial ratios found here. 

In the previous section, we described the effects on measured 
source sizes of convolving the Perseus data to lower resolution, in 
order to simulate observing it at the distance of W3. We have also 
examined the effect of this artifically reduced resolution on the ap- 
parent virial ratio for the eleven resulting convolved cores whose 
effective radii are listed in Table [9] The results are shown in Fig- 
ure [T6] where it can be seen that the effect of reduced resolution, 
and so of greater distance, is to increase the apparent virial ratio of 
sources. In this case, the aggregated Perseus region sample has a 
median virial ratio significantly closer to unity (0.67) than the na- 
tive resolution Perseus sample. The explanation of this shift appears 
to be that the increase in mass included within the larger effective 
beam is greater than the corresponding increase in velocity disper- 
sion or radius created by including additional sources. 



4 DISCUSSION 
4.1 General Results 

We find reasonable agreemen t between our results fo r ammonia 
sources in Perseus and those o f lRosolowskv et al.M2008h who com- 
pleted a single-pointing survey of ammonia observations, mainly 
toward submillimetre peaks in the Perseus region. A comparison of 
in tegrated intensity measure ments at the positions of sources listed 
in iRosolowskv et all {2008) show consistency between the differ- 
ent observations with a median difference between the values of 
15%, consistent with expected uncertainties in our observational 
techniques. Two sources (P24 and P49) show differences of greater 
than 50% at 64% and 56% respectively, although these sources are 
in the lowest 30% of our sample in terms of signal-to-noise ra- 
tio. Overall, there is a rough anti-correlation between the signal- 
to-noise ratio of a sourc e and the difference foun d between our 
observations and those of IRosolowskv et all j2008l) which is what 
one might expect for any sample subject to measurement errors. 
When comparing whole cores in our sample to values presented by 
IRosolowskv et alj d2008l) . we find also find generally consistent re- 
sults. However, there is fdisagreement in some cases, likely due to 
the fact that values presented here are integrated over the entire re- 
gion of detected emission rather than at single-pointing positions. 
There are few ammonia observations of the W3 region on the scale 
of our coverage, though there are several observations of the area 
around W3(OH) dWilson et alJI 19931 ; [Zeng et alj|l984» as well a s 
more complete coverage of W3 Main by Tieft runk et alj dl998b . 
Where our observations overlap with those available to us from the 
literature the resulting physical parameters are in good agreement. 

In terms of the NH 3 -derived source parameters, the average 
values and distributions show significant differences between the 
Perseus and W3 regions. In general terms these differences may be 
largely accounted for by the difference in distance between the two 
regions. Lower filling factors and optical depths are expected for 
regions at a greater distance due to the greater level of unresolved 
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structure within our beam. Also, if core emission is extended, this 
may account for the differences seen in observed linewidth, with 
a greater amount of material being captured by the beam in W3 
compared to Perseus, consistent with more sources being contained 
by the beam and/or a greater velocity dispersion associat ed with 
a grea ter scale of emission, as per the scaling relation of lLarsor] 
dl 98 lh . The radii covered by the GBT beam correspond to 0.04 pc 
at the distance of Perseus and 0.29 pc at the distance of W3, indi- 
cating that, while both sets of observations are capable of resolv- 
ing dense cores at the typically accepted sizes (0.1-0.5 pc), there 
may be a distinct difference in the type of structure being traced in 
Perseus as compared to W3. Convolution of the physical parame- 
ter maps derived from the Perseus data to a resolution comparable 
to that of W3 indicates that differences in distance are not able to 
explain the difference in physical parameters completely, resulting 
clump sizes are significantly smaller than native resolution obser- 
vations, indicating that there is some real difference in the dense 
structures in the two regions, with larger clusters of objects in W3 
than in Perseus. Higher average filling factors in the Perseus sam- 
ple indicate that we are resolving the structure within Perseus to a 
finer level than within W3. 

There are clear correlations between submillimetre and NH 3 
detections, although correlations in brightnesses of the two tracers 
are imperfect. The reasons for this are not clear, being reliably op- 
tically thin, submillimetre continuum emission has uncomplicated 
radiative transfer and thus is often considered a definitive tracer of 
column density. Line emission, conversely, has relatively complex 
radiative transfer and excitation processes requiring knowledge of 
or assumptions about optical depth and the relative population of 
energy levels. 

The NH 3 molecule is, like all molecules, subject to depletion 
in dense, cold gas by freezing out of the gas phase onto the sur- 
faces of dust grains. This has the effect of significantly altering 
the abundance in a somewhat unpredictable way. Such an effect 
may lead to the disparities between dust and gas emission seen in 
our observations, with the densest portions of star-forming cores 
being traced by submillimetre emission but with different chemi- 
cal specie s tracing shells o f differing chemical abundances around 
this peak JPon et alll2009h . However, nitrogen-bearing molecules 
such as NH 3 are usually predicted by models to be slow to freez e 
out, even at high densities < IChamlevll997llBergin & LangeJ 199% . 
Also, the star formation state of our observed cores makes an early- 
time depletion scenario less plausible, 83% of the NH 3 sources in 
our combined W3 and Perseus catalogue, which have an associated 
submillimetre source, are found to contain at least one 22\im source 
from the WISE catalogue JWrightetalj2010h within their 50% flux 
contour. This suggests that nearly all of the detected sources are al- 
ready star-forming and internally heated, though, any the extent of 
any heating is likely variable and may be limited. 

The value of the fractional abundance of ammonia has been 
found to have a large range in the literature, with extreme values 
ranging over four orders of magnitude. Recent val ues over large 
self-consistent samples have considerable scatter (see lWienen et alj 
|2012| ; |Punhametai]|201lh . although so far there is no well- 
established correlation of the variation of ammonia abundance with 
any physical parameter. Fractional abundance ratios found here by 
comparing integrated NH 3 and submillimetre column densities are 
in a comparitively narrow range, with a range of a factor of two 
of the mean for sources in W3 and only 30% in Perseus. These 
values are higher by between a factor of two up to an order of 
magnitude tha n those generally found in the literature at% ~ 10-7 
(c.f. x ~10~ 8 . iTafalla et al.ll2004l : iFriesen et alj|2009l) . though re- 
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cent results based o ver large survey samples dWienen et af]|2012l ; 
lDunhametal]|201lh are in better agreement. We find abundance 
ratios of ammonia to be considerably higher in W3 than in Perseus 
with a large degree of confidence, indicating that there is true vari- 
ation of this metric from region-to-region. Given this variation, as 
well as the variation found even within cores and the suggestion 
of an anti-correlation of submillimetre column density with ammo- 
nia abundance it is clear that determining NH3 abundances using 
single-point data is not advisable. 

4.2 Association of Gas and Dust Emission 

In general, submillimetre sources are found to be associated with 
NH3 emission. Submillimetre source peaks are found within NH3 
source boundaries and the detection rate of NH3 associated with 
submillimetre cores is high. However, the correlation between the 
two tracers is not overly strong, peaks i n NH3 and submilli metre 
emission do not coincide (c.f. the study o flFriesen et alj2009l of the 
Ophiuchus region) and the overall contours of emission are poorly 
correlated in a significant number of sources. In addition, an anti- 
correlation exists between peaks of the column density measured 
through submillimetre emission and the fractional abundance of 
ammonia (see Section [3~5l > with the anti-correlation stronger on av- 
erage in our W3 sample than for sources in Perseus. These results 
imply that NH 3 emission and the submillimetre continuum do not 
reliably trace the same mass component in detail, suggesting varia- 
tions in excitation conditions or deviations from LTE within cores. 
Explicitly, an inequality of dust and gas temperatures invalidate the 
assumption of local thermal equilibrium while mis-assigned tem- 
peratures create disproportionate errors in gas/dust associations. 
Biased temperature estimates affect the predicted NH 3 and submil- 
limetre column densities in opposite senses and the measured flux 
from a fixed column density is directly dependent on T in the sub- 
millimetre c ase (see Equat ion!]} and inversely in the NH 3 case (see 
Figure 3. in lMorgan et al.l j2012h ). This is likely to be particularly 
relevant when determining fractional abundances of ammonia. In 
the case that densities in the observed regi ons are lower than the 
critical density of ammonia (~4 x 10 cm 3 , lMaretetalj2009l) and 
NH 3 is sub-thermally excited then we may expect NH 3 emission 
to be fainter than predicted in a LTE analysis and so low density 
sources may be responsible for the submillimetre sources we find 
with no corresponding ammonia detection. 

Figure [T7] shows the relationship of variations found in am- 
monia fractional abundance with the submillimetre column density 
found for the corresponding sources. It can be seen that correla- 
tions exist between the two measurements, the Pearson correlation 
coefficients of the relationships are 0.73 and 0.32 for Perseus and 
W3 respectively which reveal a moderately strong correlation for 
Perseus and a marginal one for W3. The slopes of line of best 
fit on Figure [T7] show a difference in the overall relationship be- 
tween fractional abundance variation and submillimetre column 
density when compared between the Perseus and W3 regions, with 
a steeper slope found in the Perseus region. The relationship of 
fractional abundance of ammonia to submillimetre column density 
is complex as, not only may the dust opacity factor (k v ) inherent in 
the derivation of submillimetre column density vary from region- 
to-region, it may actually vary considerably within regions. This 
factor (as well as other possible influences upon the radiative trans- 
fer mechanisms of the ammonia molecule in star-forming regions) 
is dependent upon the chemistry and grain composition of the re- 
gion in question. Dust grains are expected to grow in cold, dense 
regions with potential coagulation and/or formation of ice mantles. 



This subject requires a significant level of analysis on its own terms 
and so a more thorough analysis of the relationship between frac- 
tional abundance and submillimetre column density, particularly 
with reference to internal variation within cores, will be performed 
in a subsequent publication. 

1 

4.3 Virial Ratio 

In Section [3~6l we described the determination of virial ratios for 
each region. The distribution of virial ratios in the two samples is 
similar, though with a lower median in Perseus than in W3 and 
an average close to unity in W3. The large numbers of sources 
with lower values may represent cores which have formed stars and 
started to be disrupted, or where the bound part of the core has col- 
lapsed and disappeared from the gas and dust tracers but its grav- 
itational potential energy remains. The small scatter of values to 
above unity is expected since cores with masses significantly above 
their virial mass would collapse on a free-fall time. Selection effects 
thus result in core populations with virial ratios clustering around 
unity, the lo w values of the vi rial ratio in our Perseus sample are 
explained bv lKirk et alj J2007h who showed that the low values of 
virial ratio found in Perseus may be compatible with star formation 
through the inclusion of external pressure which they modelled in 
a Bonner-Ebert sphere analysis. As star-formation activity (i.e. the 
presence of 22 micron IR sources) is observed in the majority of 
the cores in our sample, it is very likely that almost all of the cores 
are indeed bound. The difference in virial ratios in Perseus and W3 
could suggest that a significant difference exists in the star forma- 
tion environment between the regions in that external pressure is 
not particularly significant in the latter. However, the convolution 
of our observations of the Perseus region to mimic placing Perseus 
sources at the distance of W3 has the effect of shifting virial ra- 
tios closer to the values found in the W3 region, implying that the 
value of virial ratio found for a given source is a function of spatial 
resolution. Another interpretation is that star-forming clumps are 
bound on large scales but cores are pressure-confined on smaller 
scales, these possibilities are not mutually exclusive. 



5 SUMMARY AND CONCLUSIONS 

This study was designed to assess the way in which observable pa- 
rameters associated with gas and dust emission are related, a further 
goal was to compare those results between a region of high-mass 
star formation and a region of low-mass star formation. In terms of 
source brightness, our observations of ammonia in the Perseus and 
W3 star-forming regions are consistent with previous observations, 
with differences explainable by measurement error. 

There are clear differences in the measured physical parame- 
ters of dense cores in the Perseus and W3 samples. Most of these 
differences can be attributed to resolution effects introduced by the 
considerably larger distance of the W3 cloud. For example, lower 
filling factors in W3 indicate additional structure within the beam 
and larger linewidths are expected from larger scale emission. Sev- 
eral parameter differences are not so easily explained this way, 
however. Source sizes in the Per data are significantly smaller than 
in W3, even when convolved to compensate for distance, indicating 
real differences in source structure between the two regions, with 
larger clusters of objects in W3. 

A general correlation of submillimetre and ammonia emission 
has led to many studies using the two tracers in combination to 
derive, e.g. NH3 fractional abundances and virial ratios. However, 
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Figure 17. The relationship of the variation in fractional abundance of ammonia with submillimetre column density in our Perseus (red) and W3 (blue) 
samples. Lines of best fit are overplotted. 



we have shown that the correlation is imperfect in detail, with sig- 
nificant scatter in the positions of peak surface brightness and an 
anti-correlation of submillimetre-traced column density and differ- 
ential NH3 abundance within individual cores. This shows that a 
full understanding of the relationship of ammonia to the cold dust 
component traced by submillimetre emission is yet to be attained 
and that mapping of NH 3 emission is preferable to single-pointing 
measurements for determination of integrated source physical pa- 
rameters. 

This lack of detailed correlation may indicate variations in 
NH3 abundance or excitation conditions within star-forming cores. 
While ammonia is expected to be subject to depletion in late-time 
chemical models, we do not expect this effect to be widely appar- 
ent in our sample due to the large fraction of cores with associ- 
ated infrared sources, indicating internal heating. Assumptions of 
local thermodynamic equilibrium and/or the thermal coupling of 
dust and gas within may also not be reliable. 

Cases in which ammonia emission is not consistent with sub- 
millimetre emission (i.e. one tracer is detected but the other is not) 
are explicable by observational effects and/or poor signal-to-noise 
in the undetected tracer. 

The anti-correlation between NH3 abundance and 
submillimetre-traced column density within cores varies in 
significance and slope between the Perseus and W3 regions, with 
a stronger correlation and steeper slope in Perseus. Depletion 
models of ammonia do not describe a situation easily applied to 
our observations and this phenomenon requires further study. 

Virial ratios were found to be higher in the W3 sources with 
values generally below unity in the Perseus sample. However, we 
also find that, by artificially placing Perseus cores at an equivalent 
distance to W3 by degrading the resolution, the virial ratios become 



consistent between the two samples. These results may indicate a 
selection effect, but may also imply that star-forming regions may 
be gravitationally bound on large scales and pressure-confined on 
small scales. 
A ckn owledgments 

The authors acknowledge the data analysis facilities provided 
by the Starlink Project which is run by CCLRC on behalf of 
PPARC. In addition, the following Starlink packages have been 
used: Kappa, Cupid, Gaia, Convert and Coco. We would like to 
thank the helpful staff of the Green Bank Telescope in the collec- 
tion of data used in this paper related to the project GBT10C_024. 
LKM is supported by a STFC postdoctoral grant (ST/G00 1847/1) 
and DJE is supported by a STFC PhD studentship. This research 
would not have been possible without the SIMBAD astronomi- 
cal database service operated at CDS, Strasbourg, France and the 
NASA Astrophysics Data System Bibliographic Services. 



REFERENCES 

Aguirre J. E. et al„ 201 1, ApJS, 192, 4 

Allsopp J., 2012, PhD thesis, Liverpool John Moores University 

Bergin E. A., Langer W. D., 1997, ApJ, 486, 316 

Charnley S. B., 1997, MNRAS, 291, 455 

Curtis E. I., Richer J. S., 2010, MNRAS, 402, 603 

Di Francesco J., Hogerheijde M. R., Welch W. J., Bergin E. A., 

2002, AJ, 124, 2749 
Dunham M. K., Rosolowsky E., Evans, II N. J., Cyganowski C. 

Urquhart J. S., 201 1, ApJ, 741, 1 10 
Dunham M. K. et al., 2010, ApJ, 717, 1 157 



30 L. K. Morgan et al. 



Friesen R. K„ Di Francesco J., Shirley Y. L., Myers P. C, 2009, 
ApJ, 697, 1457 

Hatchell J., Fuller G. A., Richer J. S., Harries T. J., Ladd E. F, 

2007, A&A, 468, 1009 
Hatchell J., Richer J. S., Fuller G. A., Qualtrough C. J., Ladd E. F, 

Chandler C. J., 2005, A&A, 440, 151 
Ho P. T. P., Townes C. H., 1983, ARA&A, 21, 239 
Keto E., Caselli P., 2008, ApJ, 683, 238 
Kirk H., Johnstone D., Tafalla M„ 2007, ApJ, 668, 1042 
Lada C. J., Elmegreen B. G, Cong H.-I., Thaddeus P., 1978, ApJ, 

226, L39 

Landsman W. B., 1993, in Astronomical Society of the Pacific 
Conference Series, Vol. 52, Astronomical Data Analysis Soft- 
ware and Systems II, Hanisch R. J., Brissenden R. J. V., Barnes 
J., eds., p. 246 

Larson R. B., 1981, MNRAS, 194, 809 

Maret S., Faure A., Scifoni E., Wiesenfeld L., 2009, MNRAS, 
399, 425 

Moore T. J. T, Bretherton D. E., Fujiyoshi T., Ridge N. A., All- 
sopp J., Hoare M. G., Lumsden S. L., Richer J. S., 2007, MN- 
RAS, 379, 663 

Morgan L. K., Figura C. C, Urquhart J. S., Thompson M. A., 

2010, MNRAS, 408, 157 
Morgan L. K., Moore T. J. T., Eden D. J., 2012, MNRAS 
Pineda J. E., Rosolowsky E. W., Goodman A. A., 2009, ApJ, 699, 

L134 

Pon A., Plume R., Friesen R. K., Di Francesco J., Matthews B., 

Bergin E. A., 2009, ApJ, 698, 1914 
Reid M. A., Wadsley J., Petitclerc N., Sills A., 2010, ApJ, 719, 

561 

Roberts D. A., Crutcher R. M., Troland T. H., 1997, ApJ, 479, 318 
Rosolowsky E. W., Pineda J. E., Foster J. B., Borkin M. A., Kauff- 

mann J., Caselli P., Myers P. C, Goodman A. A., 2008, ApJS, 

175, 509 

Schuller F. et al., 2009, A&A, 504, 415 

Smith R. J., Clark P. C, Bonnell I. A., 2008, MNRAS, 391, 1091 
Tafalla M., Myers P. C, Caselli P., Walmsley C. M., 2004, 

Ap&SS, 292, 347 
Thompson M. A., White G. J., Morgan L. K., Miao J., Fridlund 

C. V. M., Huldtgren- White M., 2004, A&A, 414, 1017 
Tieftrunk A. R., Megeath S. T, Wilson T. L., Rayner J. T, 1998, 

A&A, 336, 991 

Ward R. L„ Wadsley J., Sills A., Petitclerc N., 2012, ApJ, 756, 
119 

Wienen M., Wyrowski F, Schuller F, Menten K. M., Walmsley 

C. M., Bronfman L., Motte F, 2012, A&A, 544, A146 
Williams J. P., de Geus E. J., Blitz L., 1994, ApJ, 428, 693 
Wilson T. L., Gaume R. A., Johnston K. J„ 1993, ApJ, 402, 230 
Wright E. L. et al., 2010, AJ, 140, 1868 

Zeng Q., Hermsen W., Wilson T. L., Batrla W„ 1984, A&A, 140, 
169 



APPENDIX A: ONLINE MATERIAL 



L. K. Morgan et al. 




L. K. Morgan et al. 




L. K. Morgan et al. 




W09 



Dust and Gas in Star-Forming Environments 

W10 



■2 



3 2 1 0-1-2-3 
Arc Minutes 

CENTER: R.A. 02 21 04.39 DEC +61 06 03.6 J2000 
W1 1 




3 2 1 0-1-2-3 
Arc Minutes 

CENTER: R.A. 02 21 52.37 DEC +61 06 22.3 J2000 
W13 



■1 



^ 



r 



3 2 1 0-1-2-3 
Arc Minutes 

CENTER: R.A. 02 22 23.83 DEC +61 06 12.2 J2000 




3 2 1 0-1-2-3 
Arc Minutes 

CENTER: R.A. 02 21 40.42 DEC +61 05 45.6 J2000 
W12 




1 -1 

Arc Minutes 

CENTER: R.A. 02 22 03.89 DEC +61 07 32.2 J2000 
W14 




3 2 1 0-1-2-3 
Arc Minutes 

CENTER: R.A. 02 23 22.01 DEC +61 12 38.5 J2000 



L. K. Morgan et al. 




40 L. K. Morgan et al. 



W31 




3 2 1 0-1-2-3 
Arc Minutes 

CENTER: R.A. 02 26 27.17 DEC +61 29 44.2 J2000 



W35/38 




ENTER: 



60 45 30 15 -15-30-45-60 
Arc Seconds 

R.A. 02 26 42.48 DEC +61 32 10.0 J2000 
W39 



-1 



3 2 1 0-1-2-3 
Arc Minutes 

CENTER: R.A. 02 26 44.35 DEC +61 29 46.0 J2000 



W33 




3 2 1 0-1-2-3 
Arc Minutes 

CENTER: R.A. 02 26 36.24 DEC +62 06 49.7 J2000 
W37 




3 2 1 0-1-2-3 
Arc Minutes 

CENTER: R.A. 02 26 44.23 DEC +61 33 24.8 J2000 
W40 



i 




t 



1*^ HJ 
in cvhwi 




3 2 1 0-1-2-3 
Arc Minutes 

CENTER: R.A. 02 26 48.38 DEC +62 08 48.8 J2000 



Dust and Gas in Star-Forming Environments 




42 L. K. Morgan et al. 




